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Abstract 

A new modelling approach that directly prescribes dynamics to the term structure of 
VIX futures is proposed in this paper. The approach is motivated by the tractability enjoyed 
by models that directly prescribe dynamics to the VIX, practices observed in interest-rate 
modelling, and the desire to develop a platform to better understand VIX option implied 
volatilities. The main contribution of the paper is the derivation of necessary conditions 
for there to be no arbitrage between the joint market of VIX and equity derivatives. The 
arbitrage conditions are analogous to the well-known HJM drift restrictions in interest-rate 
modelling. The restrictions also address a fundamental open problem related to an existing 
modelling approach, in which the dynamics of the VIX are specified directly. The paper is 
concluded with an application of the main result, which demonstrates that when modelling 
VIX futures directly, the drift and diffusion of the corresponding stochastic volatility model 
must be restricted to preclude arbitrage. 
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1 Introduction 


The Chicago Board of Options Exchange Volatility Index, which is more commonly known as 
the VIX, is an approximation of the markets’ expectation of index volatility over a 30-day time 
period. Derivatives on the VIX provide market participants with a mechanism to invest in the 
markets’ expectation of volatility, without the need for purchasing index options. Investors can 
gain exposure to volatility through the purchase of VIX futures or VIX options, which have 
become increasingly popular in recent years. In 2004, futures on the VIX began trading and 
were subsequently followed by options on the VIX in 2006. Since products are traded on both 
the underlying index and the VIX, it is desirable to employ a model that can simultaneously 
reproduce the observed characteristics of products on both indices, while remaining free from 
arbitrage. 

In this paper, a new modelling approach that directly prescribes dynamics to the term struc¬ 
ture of VIX futures is proposed. The approach is motivated by the tractability enjoyed by models 
that directly prescribe dynamics to the VIX, practices observed in interest-rate modelling and 
the desire to develop a platform to better understand VIX option implied volatilities. In the 
existing literature, many solutions have been proposed for the joint-modelling task. The ex¬ 
isting approaches can be categorised according to the assumptions made regarding the market 
of traded instruments. All existing models can generally be placed into one of the following 
categories: 

(a) The underlying index is the primary traded security. Index options, VIX futures and 
VIX options are priced relative to the underlying index. 

(b) The underlying index and a continuum of variance swaps are the primary traded securi¬ 
ties. Index options, VIX futures and VIX options are priced relative to these products. 

(c) The VIX is modelled directly. VIX futures and VIX options are priced relative to the 
VIX. 

For models that belong to category (a), the dynamics for the underlying index are specified 
under a pricing measure and the discounted price of a derivative is expressed as a local martin¬ 
gale. The square of the VIX is defined as the expected realised variance of the index and the 
discounted price of a derivative on the VIX is expressed as a local martingale under the same 
measure. 

The majority of the literature on VIX derivatives fall into this category. Zhang and Zhu 
(2006) derived an expression for VIX futures assuming Heston (1993) stochastic volatility dy¬ 
namics. Lin (2007) presented an approximation formula for VIX futures based on a convexity 
correction, which was then used to price VIX futures when the S&P500 is modelled by a Heston 
diffusion process with simultaneous jumps in the underlying index and the volatility process 
(SVJJ). A more general result was presented in Zhu and Lian (2012), who assumed the same 
dynamics for the S&P500 as Lin (2007) and derived an exact formula to price a VIX futures 
contract. The literature on VIX options is generally similar to that of VIX futures: affine 
stochastic-volatility dynamics are usually assumed for the underlying index, which enable some 
tractability to be retained when deriving option pricing formulae. A square-root stochastic 
variance model with variance jumps and time-dependent parameters was considered for the 
evolution of the S&P500 index in Sepp (2008), while option pricing formulae under the SVJJ 
dynamics were presented in Lian and Zhu (2013). Sepp (2011) and Papanicolaou and Sircar 
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(2013) attempted to capture empirically observed features of the implied skew for options on 
the VIX. The former considered a range of parametric and non-parametric models, while the 
latter employed a regime-switching stochastic-volatility model. A non-affine 3/2 plus jumps 
model was considered in Baldeaux and Badran (2014). Pricing formulae were derived for VIX 
derivatives and numerical results illustrated that the 3/2 model is a good candidate for the joint 
modelling of VIX and equity derivatives. 

For models belonging to the second approach, it is assumed that variance swap markets 
are liquid enough to justify the use of variance swaps as a model input. The extended market 
enables models to better capture term-structure features observed in both variance and volatility 
markets. Similar to the previous class of models, the VIX is defined in terms of the primary 
instruments and the discounted price of a derivative on the VIX is expressed as a local martingale 
under the same pricing measure that is used for pricing derivatives on the index. See Bergomi 
(2005), Bergomi (2008), Buehler (2006) and Gatheral (2008) for a further discussion of models 
in this class. 

A first step in proposing any model is to make assumptions regarding the class of traded 
instruments. In approaches (a) and (b), it is assumed that index instruments and possibly 
variance swaps are liquidly traded. A model is assumed for the dynamics of the liquidly traded 
instruments, which enables an expression for the VIX to be derived. By construction, these 
approaches guarantee that the assumed dynamics for the VIX are consistent with those assumed 
for the underlying securities. Hedging VIX derivatives in this framework, however, is a non¬ 
trivial exercise. Unlike other volatility related products, the VIX is not traded and it cannot be 
statically replicated, due to the non-linear transformation used in it’s definition, making hedging 
in practice a non-trivial exercise in the setup of approaches (a) and (b). 

A related modelling approach to the one proposed in this paper is that of category (c), 
in which the dynamics of the VIX are specified directly. Dynamics are assumed under a 
pricing measure and derivatives are priced as discounted expectations of their future payoffs. 
There are several examples of this approach in the literature. Grunbichler and Longstaff (1996) 
considered a mean-reverting square-root process for the evolution of the VIX and presented 
closed-form pricing formulae for VIX derivatives. Psychoyios et al. (2010) concluded that a 
mean-reverting logarithmic diffusion with jumps is supported by VIX time-series data and also 
derived closed-form formulae for VIX derivatives. A variety of model specifications were consid¬ 
ered by Kaeck and Alexander (2010). The authors evaluated the performance of a wide range of 
models for risk management and derivatives pricing applications. An empirical analysis of one¬ 
dimensional diffusions for the VIX was performed in Goard and Mazur (2013) and the authors 
concluded that a pure-diffusion 3/2 model is best suited to capture the dynamics of the VIX. 
Derivatives were then priced relative to the VIX under such a specification. Drimus and Farkas 
(2012) attempted to replicate the concept of a local-volatility surface, which originated in Dupire 
(1993), for VIX derivatives under the assumption of linear mean-reverting dynamics. The au¬ 
thors justify ignoring the dynamics of the underlying index with the claim that the VIX market 
is mature enough for the pricing and hedging of VIX options relative to VIX futures, which is 
the market practice. 

There are several complexities associated with these models. The first challenge is in ensur¬ 
ing that the joint market between the underlying index and derivatives on the VIX is free from 
arbitrage. To ensure that the markets are arbitrage free requires the derivation of restrictions 
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on the dynamics of the VIX, which is a problem that has not been addressed in the litera¬ 
ture. The derivation of restrictions that ensure no dynamic arbitrage is a well-known problem 
in other areas of finance. In interest-rate modelling, the HJM drift conditions (Heath et al. 
(1992)) ensure that there is no arbitrage when forward rates are modelled directly. Variance 
curve models for variance swap markets are analogous to forward-rate models for interest-rate 
markets. Buehler (2006) derived variance curve arbitrage conditions as well as addressing the 
problems of finite-dimensional realisations and model consistency. Many attempts have been 
made to produce similar results for option market models by directly prescribing dynamics to 
Black-Scholes implied volatilities (see for example Schonbucher (1999), Brace et al. (2001) and 
Schweizer and Wissel (2008)). The situation is much more complex for options, however, due 
to the higher dimensionality of the state space and the non-linearity of Black-Scholes implied 
volatilities. The main contributions of this paper are Theorem 5.3 and Theorem 5.4, which state 
necessary conditions for there to be no arbitrage between the joint market of VIX and equity 
derivatives. 

Another complexity associated with this modelling approach is in the appropriate specifica¬ 
tion of a market price of risk. Since the VIX is not tradable and cannot be replicated, the usual 
relationships that connect a derivative to its underlying are not typically observed. VIX futures 
are not restricted by traditional cost-of-carry relationships and VIX options violate put-call- 
parity relationships when compared to the spot. By modelling VIX futures directly, as opposed 
to the VIX itself, complexities involved with the appropriate choice of the market price of risk 
are avoided. This is similar to the comparison of short-rate models to forward-curve models in 
interest-rate modelling. 

The final motivating factor for modelling VIX futures, as opposed to the VIX directly, is the 
concept of VIX option implied volatility. To properly understand mathematical features of VIX 
option implied volatilities, a framework that connects the dynamics of the VIX futures to the 
underlying index is required. Cox and Hobson (2005) and Roper (2010) provide a comprehensive 
discussion of no-arbitrage restrictions for traditional option and implied volatility surfaces. To 
do so, the authors relied on the notion of an equivalent martingale measure. In order to extend 
these concepts to the joint market of VIX and equity derivatives, necessary conditions for the 
existence of such a measure are required. 

The remainder of this paper is set out as follows. Market conventions, the modelling frame¬ 
work and mathematical preliminaries for the paper are specified in Section 2. In Section 3, a 
general semi-martingale representation for the VIX index is derived. The dynamics are stated 
in Proposition 3.4. The representation is quite tedious and an alternate form of the diffusion 
term in a simplified setting is provided in Corollary 3.6. Section 4 is concerned with the im¬ 
plications of modelling the term structure of VIX futures directly. The dynamics of a process 
that represents the VIX as implied by the family of equations for the VIX futures is derived 
in Proposition 4.1. Section 5 is concerned with the implications of the joint modelling of the 
underlying index and the term structure of VIX futures. Restrictions on the dynamics stated 
in Section 3 and Section 4 are derived so that there is no arbitrage between the joint market of 
VIX and equity derivatives. The restrictions are stated in Theorem 5.3. The paper is concluded 
with Section 6, where an application of the main theorem is provided. The application demon¬ 
strates that by modelling the VIX futures directly, the drift and diffusion of the corresponding 
stochastic volatility model must be restricted to preclude arbitrage. 
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2 Preliminaries 


Market conventions and the modelling framework adopted throughout the paper are described 
in this section. The market definition of the YIX is first stated, followed by the theoretical 
definition of the VIX and the assumptions that are made regarding the dynamics of the traded 
instrument. 


2.1 Market Conventions 

The VIX attempts to provide investors with a model-free measure of the markets’ expectation 
of 30-day volatility of the S&P500 index. The market definition of the VIX is based on a 
representation of expected realised variance in terms of option contracts and is stated below. 


Definition 2.1. The CBOE VIX is calculated using the formula 


where 


' mrk 


VIX t (T) := 


2e r( ' T 0 ^ n, ^3 1 ( Fpcp 1 ^" 

T-tj^ph) T-t\ fc 0 1 ) 


x 100; 


( 1 ) 


T — t: is the time horizon (typically 30 days). 

r: is the risk-free interest rate that applies from time t till time T. 

ko: is the “at-the-money” strike and is given by the strike that minimizes the difference between 
put and call option prices at time t with expiry at time T. 

kji is the j th strike and kj < kj + \. 

Qj\ is the price of an out-of-the-money option at time t with strike kj and expiry at time T, 
computed as the average of the bid-ask spread. If j < 0 put options are used, if j > 0 call 
options are used and for j = 0 the average of the put and call price is used. 

I: is the set of all j. ordered by strike, for which quoted strikes exist in the market with the 
provisos that: 

— if the bid price for Qj is zero, the j (j /; 

— the summation stops if two consecutive zero bid prices are met. 

A kj'. is the central symmetric difference \{kj+ 1 — kj- 1 ), except for the first and last strikes in 
the sum where a one-sided difference, whichever of kj — kj- ± or kj+ 1 — kj is appropriate, 
is used. 


Fpcp : is the forward index level, computed using 

Fpcp := k 0 + e r( ' T ~ t \Co — Pq), 


where C*o and Pq are the mid-prices of the call and the put options at time t with strike 
ko and expiry at time T. 
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In the event that no options have exactly 30 days till expiration, as is ordinarily the case, the 
CBOE interpolates between the CBOE VIX squared calculated at the two closest maturities, 
that is, 


VIX, 


mrk 






(2) 


where t < T\ < t + 30 days < T 2 and Nt denotes the number of days from t till time T. Since 
the CBOE VIX is always defined for a 30-day time horizon, the dependence on T is omitted. 

The VIX is defined using liquidly quoted index options, however, the index itself is not 
tradable. Investors are able to take a position on the value of the VIX index via VIX futures 
and European options on the VIX, which are liquidly traded instruments. VIX futures and VIX 
options deliver a cash settlement amount that is related to the value of the VIX index at expiry. 
Traditional relationships that are usually observed between an option and its underlying do not 
hold for VIX options, as a consequence of the fact that the VIX index is not tradable. The 
relationships, however, are observed between VIX options and the VIX futures. 

For more information on the market definition of the VIX and on VIX derivatives, the 
reader is referred to the CBOE (2009) White Paper. Since the definition of the VIX is rather 
cumbersome, simplifications are adopted in the literature to enable mathematical tractability. In 
the following section, the definitions adopted in this paper and the simplifications often observed 
in the literature are presented. 


2.2 Theoretical Definition of the VIX 

In this paper, the VIX is defined in terms of European options. This is done to avoid a potential 
loss in generality, since there are no implicit model assumptions associated with such a definition. 
The definition is independent of any assumptions regarding a risk-neutral measure and absence 
of arbitrage. Upon making additional assumptions, a representation in terms of expected realised 
variance may be recovered. Until such assumptions are made, however, the representation may 
not hold. This is subtle point and it is illustrated in the forthcoming Lemma 2.5. 

Remark 2.1. The convention throughout is to use a circumflex to denote processes that are in 
terms of a fixed expiration date, while no accent is used to denote processes that are in terms 
of time-till expiry. 

Remark 2.2. For the following definition, further assumptions are required to ensure that V) is 
well behaved. For example, out-of-the-money put options need to decay fast enough as k \ 0 
to ensure that Vt does not explode. Additional assumptions are generally made to ensure that 
Vt is well defined. 

Definition 2.2. For all T > t, let 


VtiT) := 



e t (k,T;F t (T)) 

k 2 


m(dk), 


VO < t < T, 


( 3 ) 


where Ft{T) is an index future at time t expiring at time T, 0 4 (fe, T; F t (T)) is an out-of-the- 
money index option with strike k and maturity T, and m(-) is an unspecified measure. 
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Remark 2.3. The behaviour of Vt(T) as T \ t is dependent on the behaviour of ©*(•, T; F t (T)) 
as T \ i. While being of mathematical importance, this limiting behaviour is not investigated 
in this paper, since the VIX is defined for fixed T > t. 

Definition 2.3. For all i > 0 and for fixed r* > 0, the VIX is given by 

Vf.= V t (t + T*), Mt > 0. (4) 


Since t* is a fixed constant throughout, there is no loss in generality incurred by suppressing 
it as an argument to the process V. For simplicity, the dependence of the futures contract on 
t T t* is suppressed in the forthcoming analysis. 

A different choice of m(-) can be made depending on the specific application. In particu¬ 
lar, to recover a discrete-strike definition of the VIX, one would choose some form of discrete 
measure, whereas for a continuous-strike setup, one would chose the Lebesgue measure. The fol¬ 
lowing definition is motivated by the continuous-strike definition of the VIX, which is a common 
simplification observed in the literature. 

Definition 2.4. Suppose that the measure m(-) is chosen as 


m(-) = N x mi('), 

where N = 2 x 100 2 and ) denotes the Lebesgue measure. Then the VIX is given by 

and this quantity is referred to as the continuous-strike definition of the VIX. 

Under the specification given in Definition 2.4 and upon making addition assumptions re¬ 
garding the underlying index, a simpler representation for the VIX can be derived. Such a 
representation is the typical starting point for most models and is the subject of the following 
lemma. 


Lemma 2.5. Suppose that the measure m(-) is chosen as specified in Definition 2-4- Further 
suppose that there exists a measure Q, such that index futures and options are Q-martingales. 
Then the VIX can be written as 


V t = 



In 


Ft.+T* 

Ft 


T t 


( 5 ) 


Proof. For fixed t > 0, fix w £ (1. Since F t+T *(uj) € R+, Vw € D, for f(s) = ln(s) and 
X = Ft+ r * (oj), Lemma B.l yields 


ln(Ft +r .(cj)) 


ln(F 1 ) + 4(F t+T .(u,)-£) 

_ r(3±id^Vl dk 

Je k 


f E (k-F t+T *( w))+ 

Jo k 2 
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Since cj is arbitrary, the above equation is true for any oj E II. Taking conditional expectations 
with respect to Ft under the measure Q and rearranging gives 



Cf(k,t -\- t* \ Ft) dk + 


r E i 

J P(k,t + T\Ft)dk 


= -E Q [ln(F t+T .)| F t ] + In (Ft) + -{F t ~E). 


( 6 ) 


Evaluating the above equation for E = F t and multiplying both sides by 2 x 100 2 /r* completes 
the proof. □ 


The convention is to approximate the square of the VIX by the risk-neutral expectation of 
the realised variance of the underlying index. The VIX at time t is often approximated by 


V t 



x ^E«2([lnF] t+T -[lnF] t | F t ) 


r=30 days 


( 7 ) 


where [•] is used to denote quadratic variation and F t denotes an index future expiring at time 
t + T. There are many different representations of the expected realised variance of a stochastic 
processes and one of these representations, such as the logarithmic representation presented 
in Lemma 2.5, is typically chosen for the definition of the VIX. A consequence of defining 
the VIX in terms of expected realised variance is that model assumptions are implicit in the 
definition. These model assumptions are fundamental to the origins of the VIX, however, the 
market definition of the VIX at any fixed point in time is simply a function of observed market 
prices and is independent of any previously made assumptions. Lemma 2.5 may not hold, for 
example, if the index is a strict-local martingale as opposed to a true martingale. 


2.3 Index and VIX Futures Dynamics 

The modelling framework and the assumptions made regarding the class of traded instru¬ 
ments are stated in this section. The setup is rather general; the assets are assumed to be 
continuous processes that admit a finite-dimensional realisation of arbitrary size. Consider 
a continuous-time economy with trading interval [0, T*] for a fixed horizon date T* > 0. Let 
(Wt = {Wf, W f l ...., W t d )) 0 <t< T * be a (d+ l)-dimensional Brownian motion on (II, F, F, P), where 
F = (Ft)o<t<T* is the P-augmented filtration generated by W. 

Assumption 2.6. (Al) For each i = 0,..., d, let c* : [0, T*] xO^ R([_ +1 and fF : [0, T*\ xfi-> R 
be F-adapted processes, with 

T* rT* 

|/^*| ds<oo and / \a \| 2 ds < oo, P-a.s.. 

Jo 

The vector of processes (X t = (X ( t } . X f l ,..., Xf))o< t <T* € R d+1 is assumed to satisfy 


Xl = xx + 


f xW s ds+ f XI 
Jo Jo 


<y . 


d\v: 


( 8 ) 


for each i = 0, ...,d. Let F = X° , so that the underlying index future (F t = V t °)o <t<T* is 
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assumed to be a stochastic process with dynamics given by 


F t = F 0 + / F sl i° s ds+ / F s a° s -dWf. 


Jo Jo 


( 9 ) 


The motivation for such a setup is to allow for existing models, such as stochastic volatility 
models, to be examined within the proposed framework. 

Assumption 2.7. (A2) For all times t E [0, T*], the market contains call and put options 


for all strikes k E [0,oo) and time-till expiries r E (0, T* — t]. The price of an option at time 
t E [0, T*] with strike k E [0, oo) and expiry r E (0, T* — t] is derived from X* and is assumed to 


be a deterministic function, denoted by C(k,T,t,~K t ) for a call and P(k,T,t,~X t ) for a put. 

For all 0 < t < T < T*, let F l (t,T) denote the value of a VIX future at time t expiring 
at time T. The following assumption is concerned with the dynamics of the VIX futures. The 
purpose of the following assumption is to allow for a very general specification for the VIX 
futures and no interpretation or context for the measure P is initially provided. 

Assumption 2.8. (A3) For all 0 < T < T*, let v(-,T) : [0, T] x Q —>• R^ 1 and /jX (-,T) : 
[0, T] x ^ R be P-adapted processes. For all 0 < t < T < T*, the family of equations for the 
VIX futures is assumed to satisfy 





Remark 2.4. For results on the the existence and uniqueness of solutions to equations of the 
form of Equation (10), the reader is referred to Section 4.6 of Morton (1988). 

Further assume that, for any 0 < T < T*, 


T r r-T 

\/jX( u, T)| du + / |z/(it, T)| 2 du < oo, P-a.s. 


o Jo 


and that the limit 


F v (t,t) = )imF v (t,T) 


L 


•t 


L 


= F v (0,t)+ / F v (u,t) fi v (u,t) du + 


F V (u,t) v(u,t) ■ dW^, 


is well defined for all 0 < t < T*, P-a.s.. 


Given the above family of equations, one can introduce the following process that represents 
the VIX as implied by the VIX futures. 


Definition 2.9. The implied VIX is given by 


V t :=F v (Vt) 


for all t E [0, T*]. 
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Section 3 is concerned with the derivation of the dynamics of the VIX from the underlying 
index and the setup specified in (Al) and (A2) is assumed. In Section 4, the setup of (A3) is all 
that is assumed. The implied dynamics of the VIX is the focus of this section and the analysis is 
independent of any specification for the underlying index. The focus of Section 5 and Section 6 
is on the complete framework of (Al) - (A3). 


3 Deriving the Dynamics of the VIX from the Index 


In this section, a general semi-martingale representation of the VIX is derived under the mea¬ 
sure P, which is used to denote the real-world or statistical measure. A typical starting point 
for most market models is to directly specify the dynamics of the modelled quantities under 
an equivalent martingale measure. This is done to avoid the complexities involved with the 
specification of a market price of risk and with a change in measure. The reason for performing 
an analysis of the real-world dynamics is that the VIX is often used in empirical investigations, 
due to its role as an indicator for market sentiment. Starting under the real-world measure 
also avoids complexities regarding the existence of a risk-neutral measure. The ultimate goal, 
however, is to provide a framework for the pricing and hedging of derivatives, which is typically 
done under an equivalent risk-neutral measure. The existence of such a measure is not discussed 
in this section, rather, it is the subject of Theorem 5.3 and Theorem 5.4 in the forthcoming 
Section 5. 

This section is structured as follows. The Ito-Ventzel formula is first applied to derive the 
dynamics of the square of the VIX when defined in terms of options with fixed expiry. A Musiela- 
like parameterisation in terms of fixed time-till maturity is then introduced, which allows for 
the derivation of a governing stochastic differential equation for the VIX. The representation 
is presented in Proposition 3.4. In Corollary 3.5, the dynamics implied by the specification in 
Section 2.1 are stated. 

Proposition 3.1. Suppose that the setup of (Al) and (A2) in Section 2.1 are assumed, that 
is, at time t puts and calls on the index with strike k and expiry at time T are given by the 
functions P(k,T,t,~X. t ) and C(k,T,t,~K t ). Further assume that, for all T € (0, T*], 

P(;T,;-),C(;T,;-)eC and P(;T, ; ■), C(;T, ; ■) € C 1 ’ 1 ' 2 
on R + x [0, T) x R(|_ +1 , for all (k,t,x) € R+ x [0, T*) x R+ +1 , 

P(k,-,t,x), C(k,-,t,x) € C 1 


on ( t,T *], for all x € R rf+1 and 0 <t <T < T*, 


f°° 1 - 


<90 

/ -prQlk, T, t, x) m(dk) < oo, 

Jo k 

L F 

— (fc,T,t,x) 


m(dk) < oo, 


and for all x € R d+1 , 0 < t < T < T* and i,j = 0,..., d, 



<90 , 



<9 2 0 , 

L F 

fe- (w ’ x) 

m(dk) < oo, and 

i F 

(fe,T,i,x) 

OXiOXj 


m(dk ) < oo. 
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Then, for all T E (0, T*] and t E [0, T), the dynamics of the process V^(T) are given by 


„ I i r°° i 

dV?{T) = - —V t 2 (T) dt + — j ^-^-(fc, T , t, X t ) m(dk) dt 

i d poo i qA i i 

+ ~ t g l V^ T ' '■ X *» - 2(5^5?' 1 < F ' 

i d r 00 i <9 2 0 

Proof. See Appendix A. □ 

Remark 3.1. The quadratic variation term and the lack of symmetry in Equation (11) are due 
to an application of put-call parity. 

Since the VIX is defined for a fixed time horizon, it is useful to introduce a Musiela-like 
parameterisation, which is the purpose of the following lemma. One may argue that such a 
representation is unnecessary for r E (0, T* — t], due to the fact that the VIX is always defined 
for a 30-day time horizon. The representation is not entirely worthless, however, as it provides 
a starting point for the analysis of the dynamics of the VIX when the definition is based on a 
rolling portfolio (see Equation (2)). 

Lemma 3.2. For all r > 0, let 


y t (r) := V t \t + T), Vt> 0. 


( 12 ) 


Then 

1 / r°° 1 /<90 <90 \ \ 

dY t (r) = - ^ ^ + -Q^(k,T,t,Xt)\ m{dk ) - Y t (r)j dt 

i 1 d poo i an 

-^{F,F) t+ Y, l T^.r «»(*« 

* i=0 u 

1 d f°° 1 <9 2 0 

+ TrY.f WS^ (k ’ T ^ )m{ik)d( ~ X> ' Xi F 

i 1=0 u J 


P-a.s.. 
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Proof. Equation (11) and Equation (12) imply the dynamics 


Y t (r) = Yq(t) + f dY s (r) 

Jo 

ft „ ft ot>2 

= y 0 (r) + jf dVg (s + r) + J s + r) ds 

ft i ft 1 foo 1 aA 

= Yq(t)~ J —Vg(r) ds + J -J -j-j—(k,s + T,s,X s )m(dk)ds 

f ft 1 fOO 1 OQ 

+ '52J 0 tJ 0 k 2 Q~ z -{k,s + T,s,X s )m(dk)dXl 

-t d ft -I roo -) n20 

+ 2 £ l 7 i Y‘a^f k ' 5 + T ' 5 ' Xs) d <* vx ">, 

-\J Q 7Fi d{F - F) ‘ + L m\L F® (fc ’ s+T ' s ' x *) ”(*)*' 


Writing the option prices as a function of time-till expiry and differentiating under the integral, 
which is a consequence of the Dominated Convergence Theorem, yields 


Y t {r) = Y 0 (r) - 


t i roo i /O 0 

r / p[^ft T ,S)X s ) + ^ftr,s,X s )] m(dk) ds 


(r)- [ -Y s (r)ds+ j - 
Jo T Jo T 

1 Y 1 . ^ ft 1 roo 1 O0 

2 rl Ff d{F ’ F) ’ + Yf 7 1 ¥ ^.r.s.X,)m(dk) 


dXl 


1 d ft | /-CX3 1/920 

+ 2^ J 0 T l FUTTYT7. {k,r,s,X s )m(dk)d(X^X l ) ( 


i,j=0 1 


k 2 dxidxj 


Expressing the above equation in the corresponding stochastic differential equation form 
completes the proof. □ 

Definition 3.3. Let Y* = Y t (r*). Since t* is a fixed, the dependence of Y t * on t* is suppressed. 

The object of concern is the VIX, not the square of the VIX, and the following proposition 
provides the dynamics of the VIX. 

Proposition 3.4. Let v(y) := y/y, \/y > 0. Therefore V) = v(Y t *), Vt > 0, and the dynamics of 
Vt are given by 


dV t = u^Vtdt + uf ] V t d (F, F) t + 


E 4’ j) Vtd(xi, X i ) t + '£w®V t dX' t , 


i,j =0 


i =0 


(13) 
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where the drift and diffusion coefficients of V are given by the equations 


v™ = - 


( 2 ) 

u: = — 


1 l l fdQ 90 \ 

2^ + 2Fv? i F (« <*■ T< - *' x ‘> + T< - x ‘>) "**>• 


4r*V t 2 F 2 


a,j) i r 1 020 


rq = 


4r*lf J o 


k 2 dxidxj 


(k, r*,t, X t ) m(dk) — 


(i) (j) 
wf. w t 


(14) 

(15) 

(16) 


(i) 

uh = 


1 


2t*V 2 ,/q F 2 9x, 
P-a.s. ; /or all t > 0. 

Proof. By Ito’s formula, 


00 1 90 

J1 —(k,T,t,X t )m(dk), 


(17) 




dVt = w- 
dy 

where 



Therefore 



a- 

II 

to 1 

* 1 

/1 r°° i 

[v t J 0 F 2 



Id 2 - 


^ = !77( y t*) ^ + o 7172 ( y t*) d <^*>t > 


1 

2^ 


2 9y 2 

d 2 v 
dy 2 * 


4V, 3 ' 




4r*V t F 2 


4 r*V t 


i,j =0 


^-(k,T*,t,X 
00 1 <9 2 0 \ 

*<*’’**>* 


1 A / 1 90,, * v 

/ p ^-(k,T ,t,X t )m(dk) 


8 (t*) 2 V? . 4 -'_ W 0 

x y 6 z,j=0 x 



1 90 
A; 2 9a%- 


(k,T*,t,X t )m(dk)) d(X j ,X i ) 


+ 


d « c 


1 90 


2r* 17 ^ 7 0 ^ 2 <9x, 

1 i=0 u 


(k,r*,t,X t ) m(dk) dXf P-a.s.. 


Introducing it®, and w[ l \ as defined in Equations (14)-(17), completes the proof. □ 

Corollary 3.5. The dynamics of the index as specified in (Al) in Section 2.1 imply that Vt 
satisfies 


dVt 

Vt 


d 

utdi + ^w^Xiai-dW?, 
i=0 


(18) 
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where the drift and diffusion coefficients of V are given by the equations 


ut = 


E X} A - ± ^ l f (^f (*, r\t, X,) + r>, t,X, 


t) ) m(dk) - 


+ E x l x i 

i,j =0 


kfL 


1 d 2 Q 
4 t*V 2 J 0 k 2 dxidxj 


(k, r*,t, Xt)m(dfc) — 


(?) (j)' 

w\ w t 


• <d 


Ar*V t 2 


(19) 


and 


(») 

w) = 


2t*V 2 


1 90 

k 2 dxi 


(k, r*,t, Xj) m(dk). 


( 20 ) 


Proof. The result is an immediate consequence of Proposition 3.4, Equation (8) and Equa¬ 
tion (9). □ 


Remark 3.2. Given the form of the coefficients, it is not immediately obvious that there exist 
solutions to the Equation (18). Under very mild assumptions, however, V is a strictly-positive 
process and will consequently possess a semi-martingale representation with local solutions. 
To ensure that V has global solutions and a finite expectation, which are useful for practical 
purposes, additional assumptions are required. 

Remark 3.3. Many modelling approaches begin by assuming a semi-martingale representation 
for the VIX. To obtain such a representation from first principles, assumptions must be made 
regarding the option price processes to ensure that fundamental quantities exist and are well 
defined, as illustrated by the analysis performed in this section. Despite having an appearance of 
simplicity, the results of this section demonstrate that there are hidden complexities associated 
with modelling the VIX directly. In modelling the VIX directly, many implicit assumptions 
regarding fundamental quantities are made. 


3.1 An alternate representation of the diffusion term 

Proposition 3.4 provides a representation of the VIX dynamics given very general dynamics for 
the index future and index options. In what follows, an alternate representation of the diffusion 
coefficients, 


(i) i r- i d9 

given in Equation (17), is proposed. The following analysis is not affected by an application of 
Girsanov’s theorem, due to the fact that the multiplicative term is measure invariant, and is 
hence independent of any assumptions regarding market completeness and the market price of 
risk. 

The following corollary offers an alternate representation of the diffusion term and is required 
for the forthcoming analysis in Section 5. 

Corollary 3.6. Suppose that the measure m(-) is chosen as specified in Definition 2.f. Further 
suppose that there exists a measure Q, such that index futures and options are Q-martingales. 
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Then the diffusion term for the VIX can he expressed as 

^ ~ ~ ( EQ |lD(F ‘- )l Xl - *' - h '< F ‘>) ' 

Proof. The proof is almost identical to that of Lemma 2.5. Equation (6) in the proof of 
Lemma 2.5 implies that 

f°° 1 f E 1 1 

J j^C(k,T*,t,X t )dk + J ' ^P(k,T*,fX t )dk = -E®[ln(F t+T *)\F t ] + ln(F t ) + -(F t -E ). 

The conditional expectation on the right-hand side can be expressed in terms of the value of 
the process X at time t, as a consequence of the Markov property. Choosing E = Ft and 
differentiating with respect to Xj yields 


d 1 

fa-.J f2 Q ( k ’ T *,t,x)dk 


x=X, 


= (E« pn(F t+r )|X t = x] - ln(F t 


Multiplying both sides by N, dividing by 2r*Vf, and differentiating under the integral, which 
valid as a consequence of the Dominated Convergence Theorem, completes the proof. □ 


Remark 3.4. The market practice is heavily dependent on the notion of Black implied volatility, 
which is a convention that has carried across from more traditional derivative products. Tra¬ 
ditional implied volatility is well studied and it has many well-documented complexities. VIX 
implied volatilities come with many additional complexities, due to the connection between the 
underlying index and the VIX, and little progress in the way of mathematical results has been 
made. The representation of the diffusion term in Corollary 3.6 provides a potential first step 
for such an analysis. 


4 Deriving the Dynamics of the VIX from the VIX Futures 

In this section, the object of concern is the family of equations presented in Equation (10) for 
the VIX futures. The implied dynamics of the VIX is the focus of this section and the analysis 
is independent of any specification for the underlying index. To derive arbitrage restrictions on 
the joint market of VIX and equity derivatives, the dynamics of the VIX implied by the VIX 
futures must first be derived. In the following proposition, the dynamics of the VIX are derived 
from the family of equations for the VIX futures. The procedure is analogous to the recovery 
of the short-term interest rate from the forward-rate curve, which is a well known result in 
interest-rate modelling. 

Proposition 4.1. Suppose that the setup of (A3) in Section 2.1 is assumed and that the coeffi¬ 
cients pXis(t,T) and the initial VIX futures term structure, F v (0,T), are differentiable 
with respect to T, with bounded partial derivatives /j)f(t,T), vx{t,T) and F^(0,T). Then the 
dynamics of the implied VIX are given by 

Vt = V 0 + f i u V u du + f v u is(u, u) • dwl, (21) 

Jo Jo 
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where £ denotes the following process 


& = 


#(0,t) 


+ 


77 —7 [ (u,t)F v (u,t)du 

(t,t) Jo 


F v (t,t) F v (t, t) 

+ . [ (F v (u,t)v T (u,t) + Pi£(u,t)v(u,t)) ■ dWl 

F v (t,t) Jo 


Proof. Recall that the implied VIX satisfies 

rt 


V t = F v (t,t) = F v (0,t) + [ n v (u,t)F v (u,t)du+ [ F v (u,t) v(u,t) • dWl 

Jo Jo 


Therefore 


Vt = F v (0,t)+ f n v (u,t) F v (u,t) du + f F v (u,u) u(u, u) ■ dW® 

Jo Jo 


+ / (F v (u,t)u(u,t) — F v (u,u)i/(u,u)) ■ dW t 


V, 


= F 


(0, t) T f fi v (u,t) F v (u,t) du + f V u u(u,u) ■ dW^ 
Jo Jo 

+ J ^F v \u,t)u(u,t) — F v {u,u)v{u,u)'^ ■ dW^ . 

For fixed u > 0 and for each i = 0. d, 

F l ( u , t) v{u, t ) — F v (it, u ) z/(it, u) 

= J -^[F v (u,s)v(u,s)]ds 


f F v (u, s) vt{u, s) + Fiji (u, s)u(u, s) 
J U L 


ds. 


Therefore 


V) = F y (0,f) + [ n v (u,t)F v (u,t)du+ [ V u u{u,u)-dWl 
Jo Jo 


t rt r 


+ 

= F v 

+ 



0 J u L 


F v (it, s ) ut(u, s ) + Fy (it, s)n(u, s) 


V ( 


ds • dWp 


(0, i) + f /i v (u,t) F^ (u,t) du + f V u is(u,u)-dW l 
Jo Jo 


t rs r 



( 22 ) 


_F V (it, s ) vt{u, s) + Fy (it, s)v(u , s) • dW„ ds, 

/o L J 

where stochastic Fubini’s theorem has been used to interchange the order of integration. Writing 

F v {0,t) = V 0 + [ Fif (0, it) du 
Jo 

and introducing the process £ as defined in Equation (22) completes the proof. □ 
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5 Consistency Conditions 


In this section, restrictions on the dynamics of the underlying index and the family of processes 
for the YIX futures are derived. These restrictions are necessary for there to be no arbitrage 
between the joint market of VIX and equity derivatives and are referred to as consistency 
conditions. The conditions are formulated precisely in the following two definitions. 

Condition 5.1 (Cl). P(V* = V)) = 1, for all 0 < t < T*, P — a.s., where V is the process given 
in Equation (21). 

Condition 5.2 (C2). There exists an equivalent martingale measure, Q ~ P, for the underlying 
index and the VIX, such that futures and options on the index and futures on the VIX are Q- 
martingales. 

The first condition is a consequence of the restriction that the VIX and the dynamics implied 
by the VIX futures must be versions of the same process. This condition simply requires the 
two different processes that one could derive for the VIX to be in agreement. The second 
condition is a standard no arbitrage condition and it is obtained through an application of 
Girsanov’s theorem. The consistency conditions are equivalent to the forthcoming Theorem 5.3 
and Theorem 5.4. The starting point for Theorem 5.3 is an equivalent martingale measure, from 
which the drift and diffusion restrictions are derived. Theorem 5.4 starts with the drift and 
diffusion restrictions and is concerned with the existence of an equivalent martingale measure. 

Theorem 5.3. Suppose that there exists a measure Q ~ P such that futures on the index 
and futures on the VIX are Q-martingales. Then there exists a market price of risk, X, with 
/ Q T |A s | 2 ds < oo, P-a.s., such that for all T € [0,T* — r*] and t € [0, T], 

Tt = ~ At • a?, (23) 


H V (t,T) = -\fv(t,T), 


(24) 


for each j = 0,..., d, 


d / yi roo | apv 

i(M) = l va^X- T ’' t - x ‘ )m(dk)a ‘ J >’ 


with a 1 = (a 1 ’ 1 , a 1 ’ 2 ,..., a l,d ) T and 


(25) 


2 r 


1 + l X t X t w t )w t ] (<r\-°i) + X t'Yl w t x l a t 


i,j =o 


i =0 


(A t ■ v{t, t)) F v (u, t ) du — 


F v (t,t ) Jo 


ffi(0,*) 

F v (t,t ) 


F v (t 


7—7 [ (F v (u,t)u T (u,t) + Ft (u,t)v(u,t)) ■ dWl 
(t,t) Jo 


“-a.s.. 


( 26 ) 
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Proof. Let Q ~ P be an equivalent martingale measure on (fi, J-t*)- Then the existence of a 
A, with |A s | 2 ds < oo, P-a.s., is a direct consequence of the Radon-Nikodym theorem. By 
Girsanov’s theorem 


W? = W[ - / X s ds 


(27) 


defines a Brownian motion under Q. Therefore 

dF t 


= Htdt + a° • ( dW t Q + A 
= (aL° + A t ' °f) * + 


and 


d<d(k, T,t, X f ) 


no d op) I a 

TiT(*’ r - *■ x ‘) <* + E r. t, x t ) dX' t + - ]T T, t,x t )d (xKx') t 


dt 

50 


i=0 

d 


i,j =0 




50 


m (k, T, t, X t ) dt + Y, Q^(k, T, t, X t )Xf /4 dt 
2—0 

d 


i fU (E) “ ;)(E) 

+ 7 E S - S -(fcT,«,X 1 )<i<XJ,X i ) l + X; s - (k.TXWrf-dWt 


i,j =0 




*=0 


50 “ 50 

(fc, T, i, X t ) df + ^ — (A:, T, t, X t ) (ii\ + A t • at) X\ dt 


dt 


i =0 


1 ^ n2p\ ^ op) 

+ 7 E S -7 5 -(t,T,(,X t )c!<^,X i > | + ^— (k,T,t,X t )Xiai-dWi 


2 ' dxidxd 

*J=0 J 


i =0 


The no-arbitrage condition, stated as Condition 5.2 (C2), implies that, for all t G [0, T*], 

A*? = -At - at, Q - a.s., 
and that for all k > 0 and each T G (i,T*], 
rt ( op) op) 

/ — (A;,T, S ,X s )d S + ^— (k,T, 8 ,X 8 )XlfjL],ds 


2—0 


=-/„ L. E^a,».x,)xx 

*,7=0 J \ i =0 / 

) — o.s.. Introducing the change of variable T := t + r and expressing option prices as a function 
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of time-till expiry, @(-,r, •) := 0 (-,t + r, •, •), for all r £ [0, T*] and t £ [0, T* — r], gives 


j' (fc, r, a, X s ) ds + ^-(k, r, s,X s )ds + J2 f|(^ T ’ s ’ /4 ds 




3© 


i,j =0 






i=0 


( 28 ) 


for all k > 0. Corollary 3.5 and Equation (28) evaluated at r = r* imply that 


dVt_ 

V t 


1 . O-t-O-t , \^1 viyiJ'lJ) 


2r* + 4r*I4 2 




+ E 


wi 


°J ■ <4 ) + x t • E w t x 




i,j=0 


»=0 


(it 


(29) 


i=0 


The VIX futures dynamics given by the family of processes in Equation (10) satisfy, under 
the risk-neutral measure Q, 

dF y (t, T) = (p v (t, T) + A t • v{t, T)) F v {t , T) dt + F y (f, T) i/(t, T) • dW t Q . 

The no-arbitrage condition, stated as Condition 5.2 (C2), again implies that, for all t £ [0, T], 

fj, v (t,T) =-\ t ■ v(t,T), Q-a.s.. (30) 

Proposition 4.1 and Equation (30) imply that the VIX satisfies the dynamics 

dV t = itV t dt + V t i/(t, t) • dW®, (31) 


where the process £ is defined by 

F y (0,t) 1 


& = 


F y (M) F y (t 


— f ( \fi>(t,t))F v (u,t)du 
,t) Jo 


+ 


1 


(. F v (u,t) VT{u,t) + Ff ( u,t ) v(u,t)) ■ dW[ 


v ( 


F v (t, t) Jo 


Condition 5.1 (Cl) requires for the VIX and the dynamics of the implied VIX to be versions of 
the same process. Imposing the condition that the drift and diffusion coefficients in Equation (29) 
and Equation (31) must be equal completes the proof. □ 


Theorem 5.4. Suppose that pL, a, £ and v satisfy, as functions of F and F v , Equations (23)- 
(26), for all T £ [0, T* — r*] and all t € [0, T], P-a.s., for some process X, with J Q T |A s | 2 ds < oo, 


P-a.s., and 


E p 

£ ([ 

A s dWf) 

X 0 


. \J 0 

J T* 

_ 


Further suppose that there exists an adapted process F on [0, T* — r*] and a family of adapted 
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processes F v (-,T), for all T E [0, T* —t*}, on [0, T] satisfying Equation (9) and Equation (10). 
Then there exists an equivalent measure, Q ~ P, on E t *, for (F t ) 0 < t < T * and (F v (t, T)) 0 <t<T, 
for all T E [0, T*], such that futures on the index and futures on the VIX are Q-martingales. 

Remark 5.1. One such measure is given by 

<32) 

where £(•) denotes the stochastic exponential. 

Proof. The existence of a measure Q ~ P is obtained through a direct application of Girsanov’s 
Theorem. The proof of the remainder of Theorem 5.4 follows from the proof of Theorem 5.3. □ 

Remark 5.2. The martingale measure of Theorem 5.4 may not be unique due to the fact that 
there may be sources of risk that are not traded. From a practical perspective, the ability to 
replicate contingent claims is potentially of more interest than the theoretical completeness of a 
model. Theorem 5.3 and Theorem 5.4 provide an arbitrage-free representation that enables the 
direct hedging of VIX options with VIX futures contracts. 

Remark 5.3. Restrictions on both the drift and diffusion coefficients of the index and VIX 
futures are imposed in Theorem 5.3 and Theorem 5.4. The non-uniqueness of the martingale 
measure enters the theorem through the presence of A in the drift restrictions. The diffusion 
coefficient, however, is not affected by an application of Girsanov’s theorem, due to the fact that 
the multiplicative term is measure invariant. An analysis of the diffusion coefficients is hence 
independent of any assumptions regarding market completeness and the market price of risk. 

The restrictions on the diffusion terms are rather non-standard and the interpretation is 
not immediately clear, due to the presence of the term V). The forthcoming Theorem 5.6 is 
concerned with the implications of Theorem 5.3 and Theorem 5.4 on the diffusion term for the 
VIX. Before introducing the theorem, consider the following assumption regarding the diffusion 
terms. The functional form is motivated from the perspective of specifying a model for the VIX 
futures. The reason that the diffusion coefficient for the index is dependent upon V), is that this 
may be an implicit consequence from the assumed dynamics for the VIX futures. 

Assumption 5.5. The diffusion coefficient a is a function of the processes X and V, and the 
process n is a functional of the VIX futures curve F v , such that 

a t = <r(t,X t ,V t ), and v(t,T) = u{t,T, F v (t,T)). 

The partial-differential equation (PDE) in the following theorem provides a restriction on 
the choice of the diffusion coefficients, a and v, so that the joint dynamics are in agreement. 
The result is somewhat of an inverse problem; the solution to the equation is given, while the 
coefficients are unknown. 

Theorem 5.6. Suppose that the measure m(-) is chosen as specified in Definition 2.f. Then, 
for there to be no joint arbitrage in the sense of Theorem 5.3 and Theorem 5.f, for all t > 0 
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and x G R(j_ +1 , 
d 


dh 


Y x l a l ’ 3 (t,x, y/h(t,x)) — (t,x) = 2 h(t,x)is>(t,t, ^h(t,x)), j = 0,...,d, (33) 


i =0 


with h(t,Xf) = V 2 . 

Proof. Equation (25) in Theorem 5.3 states that for each j = 0, d, 

d 




i=0 


XI 


2t*V 2 Jq k 2 dxi 


7 a 7 — ( fc > r *> *> x t) rn(dk) a. 


.h3 

t I ■ 


(34) 


with cr* = (tr*’ 1 , o'*’ 2 , a i,d ) T . Recall Corollary 3.6, which states the following alternate repre¬ 
sentation of the diffusion term 

1 1 <90 N d 

- 3 —(k,T*,t,X t )m(dk)=- 


2 t*V 2 J o k 2 dxi 


2 t*V 2 dxi 


[ln(F t+T )|X t = x]-ln( J Fi) . 


Lemma 2.5 and the Markov property imply that the VIX can be represented as 


V, = -- 


-I 

T* 
d -\-1 


In 


t+T* 

Ft 


X/ = x 


=: h(t , x) 


for some function h : [0, T*] x — >• R + , and hence 




i =0 


dh 


2h{t x) ^ = ( 35 ) 


Multiplying both sides by 2h(t, x) completes the proof. 


□ 


6 Application - Proportional volatility 

In this section, an application of the main result is provided. A class of model for the term 
structure of VIX futures is considered and Theorem 5.6 is used to derive the implications of 
the model assumptions on the dynamics of the underlying index, such that the joint market 
remains free from arbitrage. The example demonstrates that there are unavoidable complexities 
involved when modelling the joint dynamics of the underlying index and VIX futures, and that 
care must be taken to avoid arbitrage. The model is a special case of models that satisfy the 
restrictions derived in Section 5, except that the risk-neutral measure Q is assumed to be fixed 
and the dynamics are directly specified under such a measure. 

The simplest non-negative class of model for the term structure of VIX futures is geometric 
Brownian motion. The influence of randomness is to shift the entire futures term structure up 
or down in a multiplicative fashion and options are priced through a straight-forward applica¬ 
tion of Black’s formula. The example provides a good illustration of the implications of the 
modelling approach, as well as a first step in an analysis of VIX implied volatilities, since these 
are calculated using Black’s formula. 
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In what follows, it is assumed that the underlying index is driven by a one-factor stochastic- 
volatility model. More precisely, let // : [0, oo) -+ R and a : [0, oo) -> K be two functions. We 
assume that the dynamics of the index process (F t ) is defined by a system of equations 


^ = y/X t dW t 

rt 

dX t = yt(X t ) dt + a(X t ) dZt , 


(36) 


where W and Z are correlated F-adapted Brownian motions under the measure Q. We assume 
that for every Xq = x > 0 and Fq = f > 0 there exists a unique strong solution (F t . X t ) to the 
system (36) such that P (X t > 0) = 1, for all t > 0, Q-a.s. In particular, 

f \p(X s )\ds + f \a (X s )\ 2 ds < oo Q — a.s. 

Jo Jo 

Corollary 6.1. Suppose that VIX futures satisfy the family of equations 


dF v (t,T) = (3(t,T)F v (t,T) dZ t , VO <t<T<T*, (37) 


for some (3 : U —>• R + where U = {( t,T ) |0 < t < T < T*}, and that the dynamics of the 
underlying index are given by the stochastic volatility model (36). 

Then, for there to be no arbitrage, fJ{t,t), p(x) and a(x) must satisfy 


f3(t,t) = 7 


and 


/i(x) = a(x) 


Ida 

2di {x) 


for some 7 € R+ and for all t > 0 and x > 0. 
Proof. By Corollary 2.5, the VIX satisfies 


V t = K -- 


N, 


In 


t+T * 

Ft 




The dynamics specified in Equation (36) further imply that 


V t = W-- 


n r f t+r * 1 r t+T * 

Ml tX’-L 


2 F? 


d (F, F ), 


Ft 


/—EQ 
! 2 T* 

- rt+T* 

J X s ds 

Ft 

/ N ™ 

■ rt+r* 


/- 

/ V s ds 

x t = 

J 2 T* 

Jt 



= : y/h(x), 

for some function h : R+ —>• R+. The function h{x) is independent of t due to the Markovian 
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structure of Xf, which can be observed by introducing a simple change of variable, 


h(x) = 





X s d-s 



X s - t ds 


X 0 = x 



f T * 


/ X r dr 

Jo 

0 * 

II 

H 

1 _ 


To proceed, introduce the function 


H(t, x ) 



r T * 


1 X r dr 

0 * 

II 

1 _ 


, 0 < t < r*, 


(38) 


such that H( 0, x) = h{x) and the function H(t , x) is the unique solution to the Cauchy problem 
(see Theorem 7.6 in Karatzas and Shreve (1991)) 


dH . . , s dH, . 1 o, ,d 2 H, . 

— (t,x) + ii(x)—(t,x) + -a 2 ( x )-^j( t , x ) + x = 0, 0 <t<r , 

H(t*,x) = 0. 


Differentiating Equation (38) and taking the limit t \ 0 implies that ¥(° + ,x) = —x. Evalu¬ 
ating the above equation at t = 0 implies that h(x) satisfies 

. .dh, . 1 9 , .d 2 h, . , . 

/i(x) — (x) + -t (7 (x) ^3 (‘ T ) = °- ( 39 ) 

Theorem 5.6 implies that <r(x) and h(x) must jointly satisfy 

f)h 

v{x)-r^{ x ) = W(ti t )h{x), (40) 


for all t > 0 and x > 0. Using the fact that h(x) is constant in time implies that f3(t, t ) = 7 , for 
some 7 > 0. Differentiating Equation (40) with respect to x yields 


da, dh. . ,^d 2 h dh . . 


and hence 


1 2 <9 2 h , 1 dh 

r ^ (x) 


9cr 

27 - & (l) 


= 7 h(x) 


9(7 

27 - & (I) 


(41) 


Equations (39)-(41) imply that 


2 q h(x) 


a(x) 


+ 7 


Ida 
2 dx 


(x) 


= 0 


(42) 
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and hence 


m 

a(x) 


+ 7 


Solving for fi(x) completes the proof. 


1 da . . 

2 & (I) = °' 


□ 


In the following corollary and example, it is shown that the assumptions of Corollary 6.1 are 
satisfied for a nontrivial class of volatilities a. 


Corollary 6 . 2 . Assume that a(x) = xao(x), where : M —> R is continuously differentiable 
and 


sup |<Jo(x)| + |<t , (:e)| < oo, x € M . 

Then, for every Xq > 0, there exists a unique global non-negative solution of the equation 


dX t = a(X t ) 




dt + a(X t ) dZ t . 


(43) 


Moreover, the process F defined in Equation (36) is positive martingale for every Fq > 0. 

Proof. The assumptions of the corollary yield the existence and uniqueness of solutions to the 
Equation (43) for every x € R. Moreover, 

X t = exp (^J ^cr 0 {X s ) Qcr' (X s ) - 7 - ^ (X s )^ ^ ds + a 0 (X s ) dZ^j . 

It is easy to check that 

E X s ds^j < 00 , 

hence (Ft) is a square-integrable martingale. □ 


Example. For a certain a > 0, consider Equation (43) with cr(x) = ay/x, for x > 0. Then the 
equation takes the form 

a 2 

dX t = —dt + a \/X t (—jdt + dZ t ). 

Recall that the equation 

a 2 r~ 

dY t = —dt + ay/Y t dZ t , Y 0 > 0, 

has unique strong solution that is instantaneously reflected at zero, hence P (T) > 0 ) = 1 for 
every t > 0. Using equivalent change of measure, Equation (43) has a unique weak solution X, 
such that P (Xt > 0) = 1 for every t > 0. Since 7 > 0 and 

E (X t ) + a 7 [ E (^X~ s )ds = E(X 0 ) + 

Jo 4 

the martingale property of (Ft) follows. 

In Corollary 6.1, the plausible dynamics for the instantaneous variance process are derived 
from the associated dynamics for the term-structure of VIX futures. A similar analysis was 
performed in Carr and Sun (2007) in the context of variance swaps. The authors propose a 
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rather general framework in which the underlying index and a single variance swap are modelled. 
Plausible risk-neutral dynamics for the instantaneous variance process are derived so that the 
dynamics of the underlying are consistent with that of the variance swap. Their restrictions are 
different to those obtained here, however, due to the fact that the object of concern is a variance 
swap, not the VIX. The authors examine the limiting case as the variance swap approaches 
expiry and, together with PDE arguments, obtain interesting restrictions from a rather general 
setup. The analogous arguments presented in Corollary 6.1 for the VIX also make use of PDE 
arguments, however, where limiting arguments were used for variance swaps, time-invariant 
arguments were used, since the VIX is defined over a fixed time horizon. 

7 Conclusion 

In this paper, a new modelling approach that prescribes dynamics to the term structure of 
VIX futures was proposed. The main contributions were Theorem 5.3 and Theorem 5.4, which 
stated necessary conditions for there to be no arbitrage between the joint market of VIX and 
equity derivatives. Not surprisingly, the restrictions are rather complex, as a consequence of 
the complexities involved in the definition of the VIX. An application of the main result was 
provided, which demonstrates that when modelling VIX futures directly, the drift and diffusion 
of the corresponding stochastic volatility model must be restricted to preclude arbitrage. This 
is similar to the well-known drift restrictions in interest-rate modelling. 

There are several directions in which the research could be extended. The analysis performed 
in this paper provides a platform for the analysis of the existing literature. For models that 
directly prescribe dynamics to the VIX, the affect of assuming a specific form for the VIX drift 
and diffusion coefficients can be assessed. Moreover, the newly developed framework provides 
a starting point for the analysis of VIX option implied volatilities, which are a fundamental 
quantity for both academics and practitioners. 

A Proof of Proposition 3.1 

Proof. For 0 < t < T < T* and fixed / > 0, let 

ff 1 „ f°° 1 „ 

G(f,T,t,X t ) = j^ j^P(k,T,t, X 4 ) m{dk) + J -^C(k,T,t,X t )m(dk) 

and write Gt(f) := G(f,T,t,X t ). Then 

V t 2 m = j^-G t (F t ). (44) 

By the assumption that Gt(-) is twice differentiable, 

SGi = P(f,T,t, X t ) _ C(f,T, t, X t ) 

Qf U) f2 /2 
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and 


d 2 G t (/) = _0_ ( P(f,T,t,X t ) _ C(f, T,t, Xf) 


dp df \ p 

Then, by Ito’s Lemma and stochastic Fubini’s theorem, 

^ 1 P 


P 


G t (f) = G 0 (f ) + 1^1 dp (k, T, s, X a ) m(dk ) + f“ ^ J q d G(k, T, s, X s ) m(dk ) 


= <?„(/) + 


o ^ Jo 

f r* i fdP 

k 2 ( <9f 



< 9 P 


/o J o 

d q 2 


(A:, T, s, X s ) ds + ^ — (fc, T, s, X s ) dX t 


i =0 


+; E pp( k ’ T ,s,x,)d{x>,x<) s I ">(<*) + 


i,j =0 
d 


dxidxj 




BC 1 1 

+ E^( t - r '»- X ') iY » + ?E | m(dk) 


i=0 


i,j =0 




G° ( /) + ^ ^ p —( k,T,s,X s )m(dk)ds + Y j J Q J o ¥ g^(k,T,s,X s )m(dk)dXi 

,T,s,X s )m{dk)d(X j ,X i ) s 



rt poo 


i=0 J 0 


t poo 


+ 


1 d 2 C 


2 fPo^ Jf k 2 dx^ 


\pU, 

i,j =0 J 


(fc, T, s, X s ) m(dife) d (X j ,X i ) s 


and 


V t 2 (T) = V 0 2 (T)- [ -J—V 2 (T)ds+[ —J dG s (F s ). 

Jo 1 ~ s Jo 1 ~~ s 

By the Ito-Ventzel formula (Lemma B.2), 


(45) 


t rF s 


G t (F t ) = G 0 (F 0 ) + [ [ 
Jo Jo 


dP(k,T, s, X s ) 


k 2 


m(dk ) 


t poo 


+ 


dC(k, T, s, X s 

P 


■ m(dk) + 


0 G S 

5/ 


(Fs) dF a 


Jf 

J 0 JF s 

+ E £ jf (^ F ‘X,s,X,) - gV„r, S ,X,)) d(X\F) t 

Pl^d^,. 
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Therefore 


G t (F t ) = Go(Fq) + 


dP(k, T, s, X .5 


- m{dk ) + 


dC(k, T, s, X s 


- m{dk) 


J (P(F„T,s,X 9 )-C(F a ,T,s,X a )') dF s 


i? y &Xj 


^ (^{Fs,T,s,X s )-^-(F s ,T,s,X s )) d(xi,F) t 


1 [* d_ ( P(f,T,s,X a )-C(f,T,s,X a 

2 Jo df[ p 


d (F, F) . 


To simplify further, observe that 


P(F t , T, t, X t ) - C t (F t ,T, t, X t ) = 0, 


E4 (Ifm.r.f.X,) X,) j <i(Xi,F) t 

j =o * V J ' J J 

P t ox 0 X0 =Ft 

= -±d{F,F) t , 


8 P(f,T,t,X t )-C(f,T,t,Xt) 
df [ P 


= - J 3 (P(/,T,t,X t )-C , (/,T,i,X t )) 
+ ^2 -§f (Hf, T, t, Xt) - C(f, T, t, Xt 

b t d J f=F t 
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for all Ft by put-call parity. Finally, 


t rFs 


G t (F t ) = G 0 (F 0 ) + 

f*t roo 
IF. 



dP(k, T, s, X s 


T 


ff 

Jo JF s 


0 JO 

dC(k, T, s, X s ) 


k 2 


-m(dk) 


k 2 


-m(dk) ff 2 d (F,F) i 


G 0 (F 0 ) + J J™ ±de(k,T,s,X a )m(dk) - ± J ^diF.F) 

ft foo i an d 

G 0 (F 0 ) + J, J 0 tf-fr( k , T , s ’Xs)m(dk)ds + Y^ 


0 r s 

d ot roo 


/ 0 J 0 

d rt roo 


i =0 Jo Jo 


i^(fc,T, S) x>w dxt 


i “ rt roo -1 £)*£) i /*£ i 

+ 2 X / 0 1 JtmF),. (46) 


Combining Equation (45) and Equation (46) completes the proof. The stochastic differential 
equation version of the result is stated in Equation (11). □ 


B Auxiliary Results 

Lemma B.l. (Baldeaux and Rutkowski (2010)) Let f : R + —>• R be twice differentiable almost 
everywhere with respect to the Lebesgue measure, so that 

rm 2 

/ \f"(k)\dk < Too, Vmj, m 2 € R+, s.t. m 2 > m\ > 0, 

J mi 

and fix y € R+. Then for any x € R, 

(,x — k) + \f"{k)\dk < Too, 

{k — x) + \f"(k)\dk < Too 

and 

f 00 ry 

f(x) = f(y) + f'(y)(x-y)+ (x — k) + f"{k)dk T / (k - x) + f"(k) dk. 

Jy Jo 

Lemma B.2. (Ito-Ventsel Formula) Let Gffx) be a family of stochastic processes, continuous 
in (t,x) € (R + x R f/ ) P-a.s. satisfying: 

(i) For each t > 0, x —>• Gffx) is C 2 from R d to R. 

(ii) For each x, {Gffx),t > 0 ) is a continuous semi-martingale 

n 

dGffx) = ^2 g\ (x)dMf 
3 = 1 

where M 3 are continuous semi-martingales, and g° (x) are stochastic processes continuous 
in (t,x), such that Vs > 0, x —>• g{{x) are C 2 maps, and Vx, gffx) are adapted processes. 
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Let X = (X 1 , ...,X d ) be a continuous semi-martingale. Then 


n pt d 

G t (X t ) = G 0 (X 0 ) + J2 9i(X s )dMi + J2 

7=1 ■'° 1=1 ' 


dG s 

dxi 


3 = 

d n ft ftJ 


(x s )dxi 


i=lj=l J[) 1 i,k=l ^ 


Proof. The version stated above is taken from Jeanblanc et al. (2009). For the original result 
see Ventzel (1965). □ 
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